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SUMMARY 


The  investigation  described  in  this  report  is  an  outgrowth 
of  a suggestion  by  D.  J.  Pastine  of  NSWC  that  the  detonation  wave 
of  a condensed  explosive  provides  a conductive  region  for  the 
deposition  of  electrical  energy  and  that  augmentation  of  detonation 
parameters  is  therefore  possible.  The  investigation  consisted  of 
preliminary  calculations  for  steady-state  detonation  in  the 
presence  of  uniform  electric  and  magnetic  fields.  Beginning 
with  the  differential  equations  governing  the  one-dimensional 
hydrodynamic  character  of  the  reacting  explosive,  we  develop 
suitable  jump  conditions  when  joule  heating  is  present  in  the 
region  of  the  reaction  zone.  It  is  apparent  from  the  results 
that  joule  heating  provides  an  energy  deposition  mechanism  that 
emulates  the  chemical  energy  release  of  the  reacting  explosive. 

In  order  to  evaluate  the  relative  merits  of  electromagnetic 
energy  deposition,  we  have  considered  an  Abel  equation  of  state 
for  TNT  as  representative  of  a condensed  explosive.  Although  the 
precise  values  of  detonation  pressure  and  velocity  depend  on  the 
equation  of  state  employed,  the  general  features  of  electromagnetic 
coupling  are  expected  to  be  unchanged.  When  the  voltage  to  which 
TNT  is  subjected  is  below  the  breakdown  voltage,  a maximum  of  12% 
increase  of  detonat*-'  -•'ressure  and  6%  increase  of  detonation 
velocity  theoretically  can  be  achieved. 

E.  Zimet  and  D.  Demske  provided  valuable  assistance  through 
discussions  of  the  validity  of  various  aspects  of  the  analysis  and 
the  accuracy  of  simplifying  assumptions.  Their  assistance  is 
greatly  appreciated.  This  work  was  supported  by  IR  funds. 

J.  W.  ENIG  ^ 
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INTRODUCTION 

The  large  energy  release  of  explosives  (typically  four 
kilojoules  per  gram  for  high  explosives)  and  the  rapidity  of  the 
detonation  process  (detonation  velocities  typically  seven 
millimeters  per  microsecond)  present  considerable  difficulty  for  real 
time  control  of  detonating  high  explosives.  In  the  case  of  gases, 
however,  the  lower  densities  and  pressures  present  a more  accessible 
regime  for  experimental  investigation.  Plett  and  Toongl 
demonstrated  the  effect  of  electromagnetic  fields  on  detonating 
gases  as  modifications  of  pressure  in  the  product  zone;  and  Mar, 
Makios,  and  Plett^  observed  wave  speeds  of  more  than  twice  the 
Chapman-Jouguet  speeds  for  hydrogen-oxygen  mixtures  when  electro- 
magnetically  driven.  Furthermore,  these  investigators  demonstrated 
a definite  electromagnetic  enhancement  of  shock  waves  in  argon. 

Thus,  the  electromagnetic  control  and  enhancement  of  condensed 
explosives  is  an  intriguing  possibility. 

Owing  to  the  overall  charge  neutrality  of  explosives,  the 
coupling  of  a detonation  process  with  electromagnetic  fields  requires 
the  microscopic  appearance  of  mobile  charge  carriers.  The  enhancement 
of  shock  waves  in  argon  shows  that  some  degree  of  ionization  must 
exist  behind  the  shock  front  even  in  a non--reactive  substance. 
Electrical  conductivities  in  the  detonation  waves  of  condensed 
explosives  were  measured  by  Cook^  in  1956;  Brish,  Tarasov,  and 
Tsukerman^  in  196n;  Hayes5  in  1965;  Ershov,  Zubkov,  and 


Iflett,  E.  G.,  and  Toong,  T.  Y.,  "Electromagnetic  Effect 
on  Flows  Behind  Gaseous  Detonations,"  AIAA  Journal,  Vol.  7, 

No.  6.  June  1969,  p.  1127.  " ' ~ ~ 

<Mar,  J.  Y.  S.,  Makios,  V. , and  Plett,  E.  G, , 
"Electromagnetic  Detonation  in  a Magnetic  Annular  Shock 
Tube,"  The  Physics  of  Fluids,  Vol.  16,  No.  12,  December  1973, 
p,  2160. 

3cook,  M.  A.,  The  Science  of  High  Explosives,  Reinhold 
Publishing  Corp.,  New  York,- 1958,  pp.  144-150. 

^Brish,  A,  A,,  Tarasov,  M.  S.,  and  Tsukerman,  V.  A., 
"Electrical  Conductivity  of  the  Explosion  Products  of  Condensed 
Explosives,"  Soviet  Physics  JETP,  Vo).  34  (10),  No.  6,  June  1960, 
p.  1095. 

5Hayes,  B. , "On  the  Electrical  Conductivity  in  Detonation 
Products,"  Fourth  Symposium  (International)  on  Detonation, 
Proceedings,  Office  of  Naval  Research,  Department  of  Navy, 
Washington,  D.C.,  October  1965,  p.  595. 


^^,^^71?F-f^' $** ''Vi5?/&£v? -J ' > VV-'^ 
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Luk 'yanchikov0  in  1974;  and  by  others.  Reported  values  of  the 
average  conductivities  vary  widely  among  investigators,  in  part 
due  to  the  extreme  difficulty  of  the  required  experimental  arrange- 
ments. But  it  is  generally  agreed  that  average  conductivities  vary 
within  the  range  of  0.1  per  ohm-cm  to  10  per  ohm-cm  for  various 
condensed  explosives  and  that  the  conduction  zone  coincides  with 
the  reaction  zone.6 


If  it  is  assumed  that  the  charge  carriers  are  electrons  and 
that  the  electrical  resistance  arises  from  collisions  of  the  electrons 
with  molecules  in  the  reaction  and  product  zoncs^  then  an  estimate 
of  the  density  of  electrons  can  be  made.  Ershov'  argues  from 
the  expression  for  the  mean  electrical  conductivity 


o »»  nee2X/mv 

19  -3 

that  the  electron  density  is  about  10  cm  for  TNT  and  about 
10 18  cm"3  for  PETN.  Here  X is  the  mean  free  path,  v the  mean 
thermal  velocity  of  the  electrons,  nQ  the  density,  and  e the  charge 
of  the  electrons. 


Such  electron  densities  cannot  be  explained  by  thermal 
ionization.  An  equilibrium  density  of  electrons  can  be  calculated 
from  the  Saha  equation, 


”■  * 2 V2  exp  (— I/kT) 

where  I is  the  ionization  potential  of  either  intermediate  or  final 
products  of  the  explosive. * 'n  order  for  the  equilibrium  densities 
to  approach  those  found  from  the  mean  conductivities,  the  ionization 
potential  must  be  about  2eV  for  TNT  and  3.5eV  for  PETN.  However, 
the  ionization  potential  for  intermediates  is  at  least  6eV  (Ershov  ). 
Thus,  the  ionization  which  accounts  for  the  conductivity  is  several 
orders  of  magnitude  above  the  expected  equilibrium  value  and  must 
be  regarded  as  originating  from  non-equilibrium  processes  within 
the  reaction  zone. 

Any  definitiv-  treatment  of  the  interaction  of  electromagnetic 
fields  with  a detonation  wave  of  a condensed  explosive  must  take 
into  account  the  mechanism  responsible  for  the  degree  of  ionization 


6Ershov,  A.  P. , Zubkov,  P.  I.,  and  Luk 'yanchikov,  L.  A., 
"Measurements  of  the  Electrical  Conductivity  Profile  in  the 
Detonation  Front  of  Solid  Explosives,"  Combustion,  Explosion,  and 
Shock  Waves,  Vol.  10,  No.  6,  November-December , 1974,  p.  776. 

f Ershov,  A.  P. , "Ionization  During  Detonation  of  Solid 
Explosives,"  Combustion,  Explosion,  and  Shock  Waves,  Vol.  11, 

No.  6,  November-December  1975,  p.  798. 

?7*Here,  ft  is  the  rationalized  Planck  constant,  approximately 
10“- ' erg-sec. 
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observed.  In  the  case  in  which  joule  heating  is  the  dominant  electro- 
magnetic contribution  (to  be  discussed  at  length  in  the  following), 
the  electrical  conductivity  of  the  detonation  wave  must  be  altered 
significantly  if  the  deposited  energy  approaches  in  magnitude  the 
energy  chemically  released.  Such  an  effort  constitutes  a significant 
research  program  in  its  own  right  and  therefore  will  be  deferred 
to  another  time.  Our  immediate  aim  is  the  characterization  of  the 
large  phenomenological  features  of  detonation  waves  of  condensed 
explosives  in  the  presence  of  electromagnetic  fields,  assuming  no 
significant  alteration  of  the  conductivity.  In  the  following,  we 
discuss  the  approximations  to  be  used,  the  energy  and  momentum  balance 
equations,  and  the  implications  of  these  considerations  to  detonation 
wave  enhancement  and  control. 


DETONATION  WAVE  CHARACTERIZATION  - FUNDAMENTAL  EQUATIONS 


The  detonation  process  of  a condensed  explosive  constitutes 
a multifluid  hydrodynamic  problem  in  which  the  various  fluid 
components  are  comprised  of  different  reacting  species.  For  the 
purposes  of  electromagnetic  coupling,  however,  it  is  sufficient 
to  regard  the  fluid  components  to  be  electrons,  ions,  and  neutral 
molecules.  Pecause  of  the  high  density  of  the  fluid  and  high 
temperature  in  the  reaction  and  product  zones,  the  collision  frequency 
is  high.  Thus,  except  for  magnetic  field  frequencies  in  excess 
of  the  electron  collision  frequency,  charge  separation  does  not 
occur.  Furthermore,  the  displacement  current  of  Maxwell's  equations 
is  assumed  negligible  in  comparison  to  other  currents  present. 

The  fluid  is  then  to  be  regarded  as  a single  conducting  fluid 
characterized  by  a (local)  conductivity  o as  well  as  the  usual 
hydrodynamic  variables  of  density,  pressure,  density,  and  temperature. 


In  addition  to  the  magnetohydrodynamic  approximation,  we  must 
determine  the  character  of  the  magnetic  field  interaction  with  the 
moving  detonation  wave.  If  the  magnetic  field  diffuses  rapidly 
through  the  conducting  fluid,  then  currents  within  the  conducting 
fluid  interact  directly  with  the  external  magnetic  field,  and  the 
effects  of  magnetic  pressure  and  induced  currents  are  negligible. 
If  x is  the  magnetic  field  diffusion  time,  V a characteristic 
velocity  of  the  fluid,  and  L a characteristic  dimension,  then  the 
magnetic  Reynolds  number 


must  be  small  for  small  magnetic  pressures.  The  diffusion  time 
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Q 

t can  be  written  in  terms  of  the  conductivity  as  (Jackson0) 

.2 


T * 


4ncL 


(2) 


Thus,  becomes 


R 


4ttoVL 


m 


,-1. 


(3) 


Using  as  characteristic  values  V ~ 7x10  cm/sec.  L ~ 10  ,cm, 
o ~ 10 2 (ohm-cm)_1(*  9xl0l3see”l)  and  c * 3xl0A0cm/sec~'L,  we  find 
R»  is  10~2.  On  this  basis,  we  will  assume  that  magnetic 
diffusion  occurs  rapidly  and  that  magnetic  pressures  are  small. 

The  equations  applicable  to  a conducting  detonating  explosive 
in  the  presence  of  an  electromagnetic  field  are  the  MHD  equations 
(Krall  and  Trivelpiece^) ; 


|f  + V-(PU)  ■ 0 , 

(4) 

9*u  i - n+  3x£ 

+ pu’Vu  * — VP  , 

(5) 

3 ■ •(*  + ¥)• 

(6) 

a - 5 ll  • 

(7) 

vxi  = 4*3  . 

c 

(8) 

Eqs.  (4)  and  (5)  are  expressions  of  continuity  and  momentum,  respec- 
tively. Eq.  (6)  is  Ohm's  law  ar.d  expresses  the  (assumed)  proportion- 
ality between  the  current  density  and  the  electric  field,  i.e., 


3 = ft 


rest* 


Eqs.  (7)  and  (8)  are  Maxwell's  equations  (in  the  cgs 


^Jackson,  J.  D.,  Classical  Electrodynamics,  John  Wiley  and 
Sons,  Inc.,  New  York,  1963,  p.  313. 

^Krall,  N.  A.,  and  Trivelpiece,  A.  W,,  Principles  of 
Plasma  Physics,  McGraw-Hill  Book  Co.,  New  York,  1973,  pp.  95-97. 
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system  o£  units)  with  the  approximation  of  negligible  displacement 
current.’’  To  these  must  be  added  an  equation  of  state. 

We  are  also  interested  in  characterizing  the  energy  of  the 
fluid.  If  e i'  the  internal  energy  per  unit  mass,  the  energy 
of  a unit  v<V>- <ne  of  fluid  is 

pe  + Jjpu2  . (9) 

Consider  the  time  rate  of  change  of  this  quantity 


f^pu2  + pe)  = |^(%pu2)  + |^{pe)  . (10) 

We  may  write  the  first  term  as 

ftiW)  ■ %u2|f  + pS-ff  . dll 

Using  the  equation  of  continuity,  we  may  write 

^ft  * “ ,*u2y*  ' (12) 


and  the  equation  of  motion  (Equation  (5))  yields 


pu* 


3u 

at 


u*VP  . 


(13) 


Eq.  (11)  can  then  be  written 


ft^PU2) 


V * (*;f>n^u) 


u*  VP  . 


(14) 


In  order  to  obtain  an  expression  for  the  second  term  of  Eq.  (10) , 
we  consider  the  thermodynamic  relation  for  a small  co-moving  element 
of  fluid 


de 


w + 


restjdt  - PdV 


(15a) 


+We  also  assume  that  dielectric  properties  have  small 
effect  on  the  detonation  process  and  thus  suppress  dielectric 
parameters  in  these  equations. 
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de  = 


¥rest) 


dt  + 


P . 
-~~dp 

P 


The  first  term  is  comprised  of  tv;o  energy  production  rate; 
the  rate  of  chemical  energy  release  per  unit  mass,  and  3*' 
the  joule  energy  deposition  per  unit  volume  of  the  moving 
of  fluid.  The  second  expression  is  valid  because  V = 1/p 
view  of  these  equations,  we  can  write 


d(pe)  = pde  + edp 

= (e  + P/p) dp  + (pw  + 34resfc)dt  . 


Thus 


|^(pe)  = (e  + P/p)f|  + pw  + 34 


rest 


By  use  of  Eq.  (4)  we  may  write  this  as 


3t 


(pe)  = - (e  + P/p)V*  (pu)  + pw  + 34 


rest 


Combining  Eqs.  (14)  and  (18)  we  have 


— (pc  + *spu2) 


- (e  -r  --)  V*  (pu) 

- V*  (Jjpu2u)  - u*VP 

+ pw  + 3'£rest  + 


This  may  be  rewritten  as 


!"£(?£  + *5P 


u2)  = - V*  £(pe  + *jpu2)uj 


+ pu*  Ve  - -V* (pu)  - u*  VP 


(15b) 


? w is 
rest  is 
element 
In 


(16) 


(17) 


(18) 


(19) 


(20) 
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From  Eq.  (15b)  we  see  that 


Ve  = L_Vp  . 


(21) 


Substituting  this  into  Eq.  (20)  we  finally  obtain 


-(pc  + ijpu2)  = - V*  £(pe  + J5pu2)uj 
- V*  (Pu)  + pw  + , 


(22) 


where  ^ refers  to  the  electric  field  in  the  laboratory  frame. 

In  order  to  appreciate  the  significance  of  the  terms  occuring 
in  Eq.  (22) , we  integrate  over  a volume  V of  the  fluid.  Then,  by 
the  divergence  theorem 


ir  I (pe  + ^pu2)  dV  = 

3t  1 v . J S 


(pe  + ijpu2)u*d§ 


- / Pu’dS  + /(pw  + 5’£)dV  . 

J S vj 


(23) 


The  left  hand  side  is  the  rate  of  change  of  energy  of  a volume  V 
of  fluid.  The  first  term  on  the  right  hand  side  of  Eq.  (23)  is 
the  energy  flux  through  the  surface  S bounding  the  volume  V.  The 
second  surface  integral  expresses  the  rate  at  which  the  surrounding 
fluid  does  work  on  the  fluid  in  V through  the  pressure.  The  volume 
integral  consists  of  volume  contributions  from  chemical  release 
(pw) , and  from  joule  heating  (J*E) . 

The  equations  which  characterize  a detonating  condensed  explosive 
in  the  presence  of  electromagnetic  fields  with  the  usual  approxima- 
tions of  low  viscosity  and  thermal  conductivity  are  taken  to  be 
the  set  Eqs,  (4-8)  and  (22)  as  well  as  an  appropriate  equation  of 
state. 


STEADY-STATE  DETONATION 

BALANCE  EQUATIONS.  Since  the  electromagnetic  contributions 
in  Eq.  (23)  appear  fn  the  same  way  as  the  chemical  release  term, 
it  is  necessary  that  these  contributions  vanish  at  some  finite  dis- 
tance behind  a detonation  front  for  steady  state  detonation  to  be 
achieved.  The  recognition  that  the  electrical  conduction  zone 
coincides  with  the  reaction  zone  (Ershov  (7))  insures  that  this 
is  the  case  since  the  current  density  J (and  thus  the  electromag- 
netic terms  in  Eq.  (23))  is  proportional  to  the  conductivity. 
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We  begin,  then,  with  a discussion  of  the  jump  conditions  for 
one-dimensional  steady  state  detonation  in  the  presence  of  an 
electromagnetic  field. 


Because  of  Eq.  (8)  we  have 

v-3  = o , 


(24) 


an  expression  compatible  with  the  assumption  that  no  charge  separa- 
tion occurs.  Since  the  conductivity  varies  within  the  detonation 
wave,  Eq-sa  (6)  and  (24)  yield  the  necessary  condition  on  the  electro- 
magnetic field, 


= 0 , 


or 


V [.(*  ♦ ¥) 

-(*  ♦¥)♦  ?•(*♦¥)= 


(25) 


(26) 


We  show  in  Fig.  1 an  idealized  one-dimensional  steady-state  detonation 
wave  as  seen  in  the  rest  frame  of  the  wave.  The  electric  and  magnetic 
fields  are  chosen  constant  and  to  have  the  directions  indicated. 

Since  the  gradient  of  a lies  solely  in  the  x direction,  we  see 
that  Eq,  (26)  is  identically  satisfied. 

The  differential  equations  (Eqs.  (4-8,  22))  are  applicable 
to  any  reacting  flow  in  the  presence  of  electromagnetic  fields 
(subject  to  the  enumerated  assumptions  we  have  made)  and,  in 
particular,  are  applicable  to  one-dimensional  flow  as  viewed  from 
a laboratory  frame,  or  as  viewed  from  a stationary  wave  frame. 

If  we  assume  that  the  chemical  energy  release  rate  can  be  character- 
ized by  a spatial  dependence  when  viewed  from  the  stationary  wave 
frame,  then  the  appropriate  equations  for  characterizing  the  wave 
are  obtained  from  Eqs.  (4-8,  22)  by  setting  all  time  derivatives 
equal  to  zero.  Eq.  (4)  becomes 


dx 


(pu  (x) ) =0  , 


or 


(27) 


pu  = const. 

We  evaluate  Eq.  (28)  at  x=0  and  find 

pu  = pcD  = ~ , 
o 

where  VQ  is  the  specific  volume  of  the  unreacted  explosive,  and 
D is  the  detonation  velocity. 


(28) 


(29) 
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Referring  to  Eq.  (6)  we  see  from  Fig.  1 that  5x§  has  a component 
only  in  the  x-direction.  Eq.  (5)  becomes 

da  . M,  . g 


From  Eq.  (29)  we  write  in  place  of  Eq.  (30) 


D_  du  _ (5x5) 


V dx 
o 


Integrating  from  x=0  to  x=°°  we  find 


n-(u  - D) 
o 


■t 


iMl_x  dx  - (Pj_  - PQ)  , 


where  u=u  (°») . Thus 


p,  + Du . p . fii  + r iMix  dx . 


1 vo  ° Vo  h c 


Also,  from  Eq.  (29)  we  have 


so  that  Eq.  (33)  becomes 


U_  _ D_ 

V ~ V ' 
1 o 


Referring  to  the  energy  equation  (Eq.  (22))  we  have 

p(e  + ^u^)u  + Pu  = pw  + 5*5  . 

In  view  of  Eq.  (29)  we  may  write  Eq.  (36)  as 


§-(e  + Hu2)  + | §-  = pw  + 54  . 
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Integrating,  we  find 


Ae  + MU2  - D2)  + - PQVo 


■H 


dx 


pw  + 


34 


(38) 


where  Ae  * e (°°)  - e (0)  . Instead  of  Eq.  (38)  we  write  the  energy 
equation  in  the  form 


Ae  + hU  + P1^1 


= HD*  + P 


dx 


pw  + 


34 


The  chemical  energy  release  per  gram  of  explosive  is 

V, 


AQ 


CH  “ D 


dx  pw  , 


(39) 


*5- 


(40) 


by  means  of  Eq.  (29) . Similarly,  we  define  an  electromagnetic  energy 
release  per  gram 


AQ 


EM 


V 

•*L 


dx  34  , 


and  write  Eq.  (39)  as 


Ae  + *U2  + P1V1  = JsD2  + PQVo  + AQch  + AQ£M  . 


(41) 


(42) 


Eqs.  (29,  35,  42)  are  the  basic  mass,  momentum,  and  energy 
balance  equations.  The  more  usual  forms  of  the  equations  can  be 
obtained  from  these.  Substituting  Eq.  (34)  into  Eq.  (35)  we  find 


u2  . 


<vo  - V 


f iMi.x  ax] 

L1  ° J 

o c 

(43) 
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Also  by  use  of  Eq.  (34) 


D2  = 


(V, 


2 r r 

fee 

[Pl  - Po  - j0 


iMi* 


dx 


o 1'  L"  JO 

Substituting  Eqs.  (43)  and  (44)  into  Eq.  (42)  yields 


♦ V,)  t & 

° 1 Jo  0 


Ae  = AQch  + AQEM  ~ *V 


+ *<P1  + V<Vo  “ Vl>  * 


x dx 


(44) 


(45) 


DETONATION  VELOCITY,  PRESSURE.  In  order  to  evaluate  the  degree 
of  influence  the  electromagnetic  field  exerts  on  a detonating 
explosive,  explicit  expressions  for  D,  P,  and  V must  be  obtained 
in  terms  of  the  initial  state  parameters  and  the  energy  release. 

The  Chapman- Jouguet  condition  and  an  equation  of  state  for  the  explosive 
and  explosive  products  are  necessary  to  complete  the  description 
and  allow  an  explicit  calculation  of  the  detonation  parameters. 

An  analysis  of  a specific  example  in  the  following  will  show  that 
strong  electromagnetic  fields  produce  relatively  small  shifts  in  the 
detonation  parameters.  Hence  an  equation  of  state  chosen  to  describe 
the  explosive  products  at  the  C-J  point  in  an  ideal  detonation  can 
be  used  to  investigate  the  degree  of  influence  of  the  electromagnetic 
field  on  the  detonation  parameters.  To  this  end  we  introduce 
(following  Eyring,  et  allO)  the  Abel  equation  of  state  for  the 
dense  product  gases 

P(V  - a)  = RT  , (46) 

where  a constant  covolume  a and  the  Chapman- Jouguet  condition 


C = V 


= u , 


(47) 


is  assumed  constant, recognizing  that  more  refined  equations  of  state 
may  be  used  without  altering  the  general  results  presented  here. 


From  Eqs.  (46,  47)  we  immediately  find 

2 2 P1 
U = V.  y 


V1  - a 


(48) 


lOEyring,  H. , Powell,  R.  E. , Duffey,  G.  H. , and  Parlin,  R. 
"The  Stability  of  Detonation,"  Reprinted  from  Chemical  Reviews, 
Vol . 45,  No.  1,  August  1949,  p.  76. 
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where  y is  the  ratio  of  specific  heats.  Using  Eq.  (48)  in  Eq.  (43) 
we  obtain 


V1  y + a8/V0 
__  - 


where  8 is  given  by 


8=1- 


12  I_ 

Pi  ‘ P1 


Y + 8 

(5x$) 


x dx  . 


0 


(49) 


(50) 


If  Cv  is  taken  to  be  the  mean  specific  heat  of  the  explosive 
in  the  range  Tp  < T < T^,  then  the  use  of  Eq.  (49)  and  the  relation- 
ship of  specif  ic"*heats  in  Eq.  (45)  yields 

(Vo  " V1)RT1 

r>  im  m \ i a i » n / rr 


where  we  have  taken 

AQ  = AQch  + AQem  - (VQ  + Vl) 

From  Eq.  (49)  we  have 

<Vo  - V t 


(3x§) 

c 


(Vx  - a)  Y ' 


vl)  , 

(51) 

dx  . 

(52) 

(53) 

and  Eq.  (52)  becomes 


r 8R 

v " 2y 


Using  the  relations 


']  T1  * °VTo  + 40  + *Vvo  - V 


Y = Cp/Cv  , Cp  - Cy  = R 


we  reduce  Eq.  (55)  to  the  form 

_2l 


T.  = 


1 Cv  8 + (2  - 8)y 


AO  + CvTq  + 


»Vo(1  - 2/V 


Y + 8 

Inserting  Eqs.  (46,  49)  into  Eq.  (44)  we  obtain 


o V 2 (y  + 8)2,.«\ 
0 = — 


V(VQ  - a) 


(54) 


(55) 


(56) 


(57) 
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which,  by  virtue  of  Eq.  (57) , becomes 

2 


D2  = 


2(V  - 1)  (Y  ± 8) 


[3  + (2  - 3) y]  (1  ~ «A0) 


AQ  + CvTq  + 


3P0V0(3.  - aA0) 
Y + 8 


. (58) 


Eqs . (49,  56,  58)  are  expressions  for  the  specific  volume  V^, 
the  temperature  of  the  reaction  products,  and  detonation  velocity 
D.  Using  the  equation  of  state  we  can  write  as 


P * ill+J), , ll_=.  1) 

*1  (Vc  - a)  [3  + (2  - 3)y 


AQ  + CvTq  + 


3P0voa  ~ «A0) 

Y + 8 


. (59) 


These  equations  characterize  the  one-dimensional  steady  detonation 
wave  in  terms  of  the  unreacted  state  of  the  explosive,  the  energy 
release  per  gram  cf  explosive,  and  the  electric  and  magnetic  field 
strengths  present. 


STEADY-STATE  DETONATION  (5=  0) . The  magnetic  field  appears 
in  the  expressions  for  the  ’final  state  parameters  and  the  the  deto- 
nation velocity  through  two  terms,  AQ  and  3.  In  the  next  section 
we  show  that  for  fields  substantially  less  than  a megagauss,  the 
approximation  that  B = 0 is  valid. 


Referring  to  Eq.  (50)  we  see  that  in  the  absence  of  5 


3=1-  PQ/P1  . (60) 

Even  without  the  added  energy  release  due  to  the  electric  field, 
the  final  pressure  of  reaction  products  is  considerably  larger  than 
the  initial  pressure  of  a condensed  explosive.  Thus,  in  Eqs.  (49, 
57,  58,  60)  we  set  PQ  = 0 and  3=1  (as  well  as  B = 0) . We 
have  then 


V,  1 + cxA, 


Y + 1 ' 


T,  = 


D 


JLl 


Cv(l  + Y) 


= nr  - i? 

(1  - aA0) 


AQ($  = 0)  + CVTQ 


AQ($  = 0)  + CvTq 


= It X - 1) 
(Vc  - a) 


AQ = 0)  + CvTq 


(61) 

(62) 

(63) 

(64) 
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Eqs.  (41,  53)  yield 


AQ($  = 0)  = Q, 


V /*° 

•■T3- 


£ dx  . 


Using  Eq.  (65)  in  Eq.  (63)  we  write 


(1  - a/V 


^ AQCH  + CvToj 


2 (Y  - 1)V 


(1  - a/V 


— j-  r 5* 

o> H 


$ dx  . (66) 


The  first  term  on  the  right  hand  side  of  this  equation  i3  the 
square  of  the  detonation  velocity  in  the  absence  of  the  electric 
field,  i.e.  D02.  Dividing  both  sides  of  Eq.  (66)  by  D02  we 


find 


2^_  = i + Za  i ^ dx 

Do2  ' D <4«ch  + CvTo>  ' 


^ r . 

»L  3-£d* 


Then 


(4QCH  + W • 


n2  D 

P_  = i + r— 
2 x T 1 D 
D 

o 


- r = o . 


The  detonation  velocity  is  thus  determined  by  a cubic  equation 
in  which  the  adjustable  parameter  is  the  joule  deposition  term  as 
given  by  T.  The  quantity  T is  a measure  of  the  ratio  of 
electrical  power  dissipation  to  a chemical  power  release  of  the 
explosive  in  the  absence  of  electromagnetic  interaction.  Once  D/D0 
is  determined  by  the  specification  of  r,  the  pressure  and  temperature 
of  the  reaction  zone  (relative  to  their  values  in  the  absence  of 
electromagnetic  interaction)  can  be  determined.  From  Eqs.  (62,  64) 
we  have 
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T./T , = 1 + Dr/D  , 

O 


P ^ = 1 + DQr/D  , 

x O '"' 


(70) 

(71) 


where  the  subscript  "oB  has  been  used  to  designate  values  in  the 
absence  of  electromagnetic  interaction. 

Approximate  Behavior.  The  effect  of  both  large  and  small 
electrical  energy  deposition  on  the  detonation  parameters  can  be 
estimated  from  the  preceding  equations.  In  the  case  of  small  T 
we  may  assume 


D/Dq  ^ 1 + 6 , <«1  . 
Then  from  Eq.  (69)  we  find 

6 = T/2  . 


(72) 


(73) 


We  may  write  for  small  r 


D/Dq  = 1 + T/2  , (74) 

Pl/Pl  3 T]/Tl  = 1 + r • (75) 

o ' o 


In  the  limit  of  large  T it  is  immediately  clear  from  hq.  (69) 
that  D/Dq  is  large  compared  to  unity.  We  can  thus  write 

^ T1/3  , T>>1  , (76) 


and  from  Eqs.  (70,  71) 


P1  *^1  2/3 

—•  = * r^/J  , r>>i  . (77) 

xO  xO 

Figures  2 and  3 provide  curves  of  the  detonation  parameters 
for  intermediate  values  of  T. 

Finally,  as  a specific  example  we  apply  the  above  equations 
to  TNT.  Table  1 is  taken  from  Eyring,  et  al^O.  in  order  to  simulate 
the  one-dimensional  character  of  the  detonation  process,  we  imagine 
a slab  of  explosive  positioned  between  two  parallel  conducting 
sheets,  A detonation  wave  is  initiated  so  that  it  propagates  in  a 
direction  parallel  to  the  sheets,  and  thus  the  wave  remains 
perpendicular  to  the  conducting  walls.  If  the  walls  are  maintained 
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at  a potential  difference  V,  a constant  electric  field  is  produced 
which  remains  parallel  to  the  detonation  wave  surface.  The 
geometrical  requirements  of  Pig.  1 are  thus  met. 

Because  the  unreacted  explosive  is  in  contact  with  the 
conducting  walls,  (in  order  to  provide  a conducting  path  for  the 
moving  detonation  wave)  arbitrarily  large  potential  differences 
cannot  be  suppc  ted  without  dielectric  breakdown  of  the  unreacted 
explosive  and  subsequent  initiation  at  points  in  front  of  the 
detonation  wave.  The  breakdown  voltage  of  TNT  has  been  found  to 
be11  3.7  x 104  volts  cm*"1.  Values  of  breakdown  voltages  for  various 
other  explosives  have  been  obtained  by  J.  W.  Forbes  at  NSWC/White  Oak 
and  are  given  in  Table  2. 

Table  1 

Detonation  Characteristics  for  TNT 
(Source;  Eyring,  et  al.,  Ref.  10) 


VD  - 0.6360  cm^gnT*  , 
Cv  = 0.326  cal  deg  gm 
y = 1.240  . , 

D0  = 6. 85x10°  cm  sec 
AQ  = 1080  cal  gm-1 
a = 0.4225  cmdgm**1 


Table  2 

Voltage  Breakdown  Levels  of  Secondary  Explosives 


Densi ty 

Distance  Between 

Breakdown 

Sample 

(qm  cm_3) 

Electrodes  (cm) 

Voltaqe  (V) 

Laboratory  air 

UT» 

.165 

5,700  ± 160 

Silicone  potting 

compound  (RTV-61) 

1.02 

.08 

20,000  ± 200 

HBX-1  (cast) 

1.69 

.20 

410  ± 55 

Type  S-l  (cast) 

1.69 

.20 

1,550  ± 90 

Comp  B (cast) 

1.62 

.20 

10,200  ± 420 

TATB 

1.84 

.20 

11,500  ± 890 

PBX-9404 

1.81 

.20 

> 24,000 

11Fedoroff,  B.  T. , 

and  Sheffield,  0.  E. , Encyclopedia  of 

Explosives  and  Related 

Items,  PA.TR 

2700,  Vol.  5,  Picatinny 

Arsenal,  Dover,  N.J.,  1972,  p.  D 1221. 
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Table  3 

Electrical  Conductivity  of  Various  Explosives 
(Source:  Brish,  Tarasov,  and  Tsukerman,  Ref.  4) 


Explosive 

Alloy  of  50% 
& 50%  RDX 


Density 
(qm  cm~3) 


Detonation 

Speed 

(xlO  cm^sec"^) 


Average 

Conductivity 

(ohnrlcni-!) 


TNT 


i 

(molded) 

1.64 

7.7 

5.0 

i 

50%  TNT  and  50% 

f 

RDX  (powdered) 

0.8 

5.7 

1.0 

RDX  (powdered) 

0.8 

6.2 

0.23 

i 

TNT  (powdered) 

0.8 

4.5 

4.0 

1 

PETN  (powdered) 

0.8 

5.2 

0.13 

>■' 

Tetryl  (powdered) 

0.8 

5.0 

0.16 

1). 

Lead  azide  (powdered) 

^0.9 

2.7 

6.2 

We  take  the  breakdown  voltage  of  TNT  as  a practical  limitation  on 
the  maximum  field  to  which  the  explosive  is  subjected  as  an  illustra- 
tive example,  bearing  in  mind  that  other  procedures  exist  for 
achieving  greater  field  strengths  within  the  conduction  zone. 

. 4 

Table  3 is  taken  from  Brish,  Tarasov,  and  Tsukerman  as  reliable 
conductivity  values  for  secondary  explosives.  We  see  that  the 
conductivity  varies  from  4 ohm~^cm~l  for  powdered  TNT  to  5 ohm"^cm~l 
for  Composition  B,  We  take  a value  of  5 ohm“lcm~l  as  representative 
of  TNT  in  the  above  equations. 


Collecting  the  results  of  the  previous  paragraphs,  we  have 

o=5  ohm  ^cm  * , , 

E a 3.7  104  volts  cm-  , (78) 


L = 0.1  cm 


where  the  conduction  zone  length  is  approximately  the  reaction  zone 
thickness  (Ershov7).  Referring  to  Eq.  (68)  we  can  write  T in 
terms  of  the  mean  conductivity  as 


° r“ir2 
r LoE 

o 


<AQCH  + cvV 


(79) 
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We  find  that  for  the  values  appropriate  to  TNT  th*t 


T = 0.13 

(80) 

and  a value  of  J 

5 -2 

J = 1.9x10  amp  cm 

(81) 

Reference  to  Fig. 

(3)  gives 

P /P.  = T /T  * 1.12  , 

x xo  1 Jo 

(82) 

and 

D/Dq  * 1.06  . 

(83) 

Thus  the  pressure  and  temperature  in  the  product  zone  are 
increased  by  approximately  12%  by  the  energy  deposition  provided 
by  the  electric  field,  and  the  detonation  velocity  is  increased 
by  approximately  6%  over  what  is  obtained  in  the  absence  of  electro- 
magnetic interaction. 

STEADY-STATE  DETONATION  (B  / 0) . An  examination  of  Eqs.  (57, 
59,  60)  as  well  as  Eg.  (50)  shows  that  the  magnetic  field  influence 
arises  primarily  through  the  ratio 


We  can  estimate  the  magnetic  field  strength  necessary  to  produce 
significant  alteration  of  the  detonation  parameters  for  the  case 
of  TNT  discussed  in  the  previous  section.  Let  Pi  be  the  product 
zone  pressure  in  the  absence  of  interaction.  Then  a magnetic  field 
necessary  to  make  the  pressure  ratio  to  be  of  the  order  of  unity 
is  given  by 


or 


B vT 


9 


(84) 


where  J is  the  mean  current  density  for  which  we  assume  the  value 
in  Eg.  (91) . Referring  to  Table  1 we  find  B iq8  gauss;  it 
appears  that  only  for  strong  magnetic  fields  can  the  detonation 
parameters  be  substantially  altered. 
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m 


1 


A more  severe  limitation  on  the  magnetic  field  strength  in 
the  applicability  of  the  preceding  analysis  is  found  on  closer 
scrutiny  of  the  mechanism  of  electrical  conductivity.  The 
conductivity  is  determined  by  the  number  density  of  the  electrons, 
their  mean  velocity,  and  the  collision  cross-section  with  ions. 

The  ions  are  themselves  screened  by  the  mobile  electrons  providing 
an  effective  cutoff  to  the  Coulomb  scattering  at  the  Debye  length 


5d  ■ 


kT. 


4ire^n 


(85) 


where  Te  is  a representative  electron  temperature  and  n§  is 
the  number  density  of  the  electrons.  In  the  presence  of  a magnetic 
field,  the  collision  cross-section  of  the  electrons  and  ions  may 
be  assumed  unchanged  only  if  the  radius  of  curvature  of  the  electrons 
in  the  magnetic  field  exceeds  the  Debye  length.  Free  electrons 
in  a magnetic  field  have  an  orbital  radius 


r 


R = 


3kT  me" 
e 

? 2 
e B 


(86) 


Equating  R and  we  find  that  B has  a value 


B J' 
max 


12TTnmc" 


(87) 
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If  we  use  a representative  value  for  n of  10  cm  , we  find 

Bmax  ^ 5x-^6  gauss  . This  is  two  orders  of  magnitude  lower 

than  the  estimate  determined  from  pressure  considerations,  even 
for  the  large  value  of  current  assumed.  We  therefore  treat  the 
magnetic  field  effects  as  perturbations  on  the  jump  conditions. 


Let  A be  given  by 


A = 


-iMl-x  a*  . 


(88) 


A is  a measure  of  the  magnetic  interaction  relative  to  the  detonation 
pressure  in  the  absence  of  electromagnetic  interaction.  By  the 
arguments  of  this  section  it  can  have  a value  as  large  as  0.1. 

Instead  of  Eq.  (50)  we  write 


8 - 1 - A 


(89) 
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valid  to  first  order  in  A (Pq/Pj  is  small  compared  to  A) . 

Each  of  the  expressions  for  the  detonation  parameters  can  be 
written  in  terms  of  their  values  in  the  ebsence  of  electromagnetic 
interaction.  From  Eqs.  (57,  59,  60)  we  write 


where 


R 


!i  = L LUEl-.- R 

T1  \2y  + (1  - y)3]  K ' 
o 


1L.  « .. . h,..±.  g ) , r 

p pY  + (i  - y)b]  k ' 
o 

°L  = fr...±.6li R 

D 2 (Y  + 1)  pY  + U - Y)g]  ' 


AQ  + CvTq 
AQCH  + CvTo  ' 

1 + IT  r (1  - d/V~)  (Y  + 3)  [2y  + 6(1  + *S*o'] 


(90) 

(91) 

(92) 

(93a) 

(93b) 


The  expression  Eq.  (93b)  is  found  f’  om  Eq.  (53)  and  from  Eq.  (49) , 
which  we  reproduce  here  for  convenience 


Y + ct3/V0 
Y + 6 


To  first  order  in  A we  find 


(49) 


x 


+ xA) 

2 (Y 
. * 


- 1)  (2Y  +(1  + a/V 
7y  + 1)  (1  - a7v^T 


(94a) 


, (94b) 


(95a) 
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y = - 


+ r 


o |A=0 


Y + 1 ' 


-1 


xr  + 2 


(y  -1)  (2y  + (1  + a/vc)) 


o JA=0 


(Y  + 1)  d - a/Vo) 


(95b) 


T, 


(1  + zA)  , 


A=0 


(96a) 


z = y + 


Y + 1 ' 


(96b) 


V, 


o JA=0 


1 + A| 


_1 _ 

Y + 1 (Y  + a/'V  ) 


(97) 


The  sign  of  A can  be  determined  by  the  choice  of  the  direction 
of  the  magnetic  field  (Eq.  (88)).  It  is  apparent  from  these  equations 
that  the  magnetic  field  can  have  an  accelerating  or  retarding  effect 
on  the  detonation  wave.  The  quantities  x and  y are  intrinsically 
negative  leading  to  an  enhancement  of  detonation  velocity  and 
pressure  by  the  magnetic  field  if  A is  chosen  to  be  negative, 
and  a retarding  effect  otherwise. 

Returning  to  the  example  of  TNT,  we  determine  the  extent  of 
magnetic  interaction  for  the  largest  current  which  can  be  established 
subject  to  breakdown  limitations,  namely  a current  density  of 
1.9x10^  amp  cm“2.  Using  the  parameters  given  in  Table  1 and 
Figure  3 we  obtain  the  values  x = -1.58,  y = -2.70,  z = -2.26. 

If,  finally,  we  choose  a value  of  A of  -0,1  corresponding  to  a 
magnetic  field  of  5x10^  gauss  (with  an  appropriate  direction 
chosen) , we  see  that  the  detonation  velocity  is  further  increased 
by  approximately  16%  and  that  the  detonation  pressure  is  further 
increased  by  approximately  27%  above  the  values  obtained  by  joule 
heating  alone.  Note  that  both  the  detonation  velocity  and  pressure 
changes  occur  in  the  same  direction  for  a given  A (x  and  y are 
both  negative) . 


DISCUSSION 


We  have  obtained  the  jump  equations  for  steady-state  detonation 
of  condensed  explosives  in  the  presence  of  stationary,  uniform, 
electric,  and  magnetic  fields.  Final  state  pressures  and  temperatures 
and  the  detonation  velocity  have  been  characterized  for  the  model 
using  an  Abel  equation  of  state.  The  prescription  can  be  used  with 
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more  realistic  equations  of  state,  but  it  is  expected  that  the 
qualitative  features  of  this  report  would  be  unchanged. 

We  can  conclude  from  the  specific  example  of  TNT  that  the 
control  of  the  detonation  wave  of  a condensed  explosive  by  means 
of  electromagnetic  fields  is  a real  possibility  requiring  electric 
field  strengths  on  the  order  of  104  volts  cm~l  and  magnetic 
field  strengths  on  the  order  of  106  gauss. 

The  analysis  presented  here  must  be  regarded  as  preliminary 
because  of  the  approximations  used.  A more  substantive  analysis 
should  treat  the  multifluid  character  of  the  reacting  explosive; 
the  treatment  of  the  reaction  and  product  zones  as  a three-component 
fluid  of  electrons,  ions,  and  neutral  particles  would  provide  a 
more  realistic  f ramewoi x for  understanding  the  details  of  electro- 
magnetic interaction  with  a detonating  explosive.  The  expression 
of  Ohm's  law  used  in  this  Report  is  at  best  a crude  approximation 
for  this  problem  and  can  be  supplanted  by  a more  accurate  expression 
derived  from  the  species  equations  of  motion. 

The  conductivity  of  the  detonating  explosive  has  been  assumed 
to  be  unchanged  by  the  presence  of  the  electromagnetic  field.  This 
was  assured  in  part  by  the  restriction  on  the  allowable  magnetic 
field  strength.  In  view  of  the  fact  that  the  current  is  induced 
by  the  local  fields,  the  joule  heating  depends  quadratically  on 
the  electric  field  strength,  and  the  magnetic  terms  are  proportional 
to  (with  proportionality  factor  o)  and  depend  bilinearly  on  E 
and  B.  Thus,  deviations  of  the  conductivity  should  only  arise  in 
higher  order  considerations. 

One  other  assumption  has  been  made  with  regard  to  the  dielectric 
properties  of  the  explosive.  It  has  been  assumed  that  no  induced 
polarization  occurs  either  in  the  unreacted  explosive  or  in  the 
reaction  zone.  Since  the  stored  energies  in  induced  polarization 
are  small,  this  does  not  change  the  jump  conditions,  but  it  does 
affect  the  field  distribution  within  the  reaction  zone  if  the  field 
is  generated  by  an  external  capacitor  plate  configuration.  This 
is  an  added  feature  which  can  be  incorporated  in  the  multifluid 
treatment. 

Finally,  it  is  worth  mentioning  that  the  maximum  field  strength 
to  which  the  detonation  wave  could  be  subjected  is  limited  by  the 
breakdown  voltage  of  the  unreacted  explosive.  A similar  limitation 
could  exist  for  the  product  zone.  Korol'kov,  Mel'nikov,  and 
Tsyplenko  12  have  shown  that  the  breakdown  voltages  of  the 
detonation  products  of  PETN  and  RDX  are  of  the  same  order  as  the 
unreacted  explosives  themselves.  Hence,  the  breakdown 
voltage  limitations  are  unchanged. 

II 2 Korol'kov,  V.  L. , Mel'nikov,  M.  A.,  Tsyplenko,  A.  P. , 

"Dielectric  Breakdown  of  Detonation  Products,"  Soviet  Journal 
of  Technical  Physics  (JETP) , Vol.  19,  No.  12,  June  1975, 
p.  1569. 
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